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i. ■ Abstract 

^! Deep saline aquifers are promising geological reservoirs for CO2 seques- 

t/5 I tration, if they do not leak. The absence of leakage is provided by the caprock 

integrity. However, CO2 injection operations may change the geomechanical 

^- stresses and cause fracturing of the caprock. We present a model for the 

^ . propagation of a fracture in the caprock driven by the outflow of fluid from a 

, P <! low-permeability aquifer. We show that to describe the fracture propagation, 

one needs to solve the pressure diffusion problem in the aquifer and in the 

fracture. The difference between permeability of the aquifer and permeabil- 



cn ■ ity of the fracture allows us to solve the problem analytically, by virtue of its 



self-similarity. We obtain an analytical expression for the fracture length as 
■^ ! a function of time. Using the obtained expression we estimate the length of a 

^T) \ hypothetical fracture at the In Salah CO2 injection site to be of the order of 

O ■ 25 — 35 m within 10 years after fracture initiation. This result shows that if 

a fracture appears it is likely that it will become a pathway for CO2 leakage. 

rS • 1. Introduction 



Use of fossil fuels for satisfaction of current energy needs has an inherent 
waste product - carbon dioxide. Since the beginning of the technological 
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revolution the amount of CO2 released in the atmosphere has grown mono- 
tonically, causing a substantial increase of its concentration. Within the last 
decade, a significant effort has been expended on identifying ways to avoid 
CO2 release in the atmosphere utilizing CO2 sequestration. Various geo- 
logical formations are considered as options for long-term storage of CO2: 
depleted oil aquifers, unmineable coal seams, deep saline aquifers, etc. The 
latter are especially promising because they are widespread and have high 
capacity. 

Deep aquifers are separated from the shallow freshwater aquifers by caprock 
- a formation with extremely low permeability (often shale). When CO2 is 
injected into an aquifer, the integrity of the caprock prevents CO2 leakage. 
However, buildup of the fluid pressure caused by CO2 injection changes the 
stresses in the caprock, and can lead to reactivation of preexisting faults [l| 
or even fracturing of the caprock [2] . 

Recent studies have shown that when CO2 is injected at a temperature 
lower than the ambient temperature of the formation, additional thermal 
stresses develop around the injection well and the risk of fracturing increases 
|3|, 1^, so that even the caprock can be fractured |5|, p]. However, fracturing 
does not necessarily lead to leakage; CO2 will leak out of the aquifer only if 
the fractures are long enough to reach an abandoned well |7| or connect to a 
network of natural fractures [Sy . The initial length of a fracture can be small, 
but under high pressure the fracture may propagate. Therefore, the rate of 
fracture propagation and the characteristic length of fractures are crucial for 
assessing the possibility of CO2 leakage from a deep aquifer. 

Multiple physical processes are involved in fluid-driven fracture propaga- 
tion: fracture mechanics, flow in the fracture and flow in the porous aquifer. 
However, when an aquifer has low permeability, the fluid outflow from it is 
slow and therefore it is the rate-limiting process for fracture propagation. 
This feature allows us to reduce the problem of fracture propagation to solv- 
ing the pressure equation in the aquifer and inside the fracture. 

Even when considering the pressure diffusion only in the vertical direc- 
tion, the problem is non-trivial since it includes two regions with different 
permeabilities and has a moving boundary condition at the fracture tip, due 
to the fracture propagation. However, we show here that when the permeabil- 
ity of the aquifer is significantly lower than the permeability of the fracture, 
the solution to the diffusion problem is self-similar. In this case, we reduce 
the original partial differential equation to an ordinary differential equation, 
which we solve analytically. Our model allows us to obtain an analytical 



expression for the fracture length as a function of time. 

Using our analytical solution we make estimates based on the parameters 
for the Krechba aquifer (In Salah, Algeria) from |9|, |5|. This site is of 
significant technological interest because is has been used as a pilot project 
for CO2 injection since 2004. We show that although the fracture propagation 
is slow, a model fracture still extends up to ~ 25 — 35 meters within 10 years 
after initiation. On such a length scales a fracture may easily reach a leaky 
fault, a system of natural fractures or an abandoned well and become a 
pathway for CO2 leakage from the aquifer into potable aquifers or even into 
the atmosphere. 

2. Problem Formulation and Model 

We consider the physical system to consist of a porous aquifer filled with 
fluid (brine and injected supercritical carbon dioxide) and the caprock (shale) 
that constrains the aquifer from above. We assume that the aquifer has 
relatively low permeability (~ 10 — 100 mD), which is the case, for example, 
for the sandstone aquifer at the Krechba field (In Salah, Algeria). Injection 
of CO2 leads to a pressure buildup in the aquifer and therefore to tensile 
stresses in the caprock. When tensile stresses exceed the caprock strength, 
the caprock fractures. We consider a single vertical 2D fracture originating 
at the boundary between the aquifer and the caprock, and we assume that 



the fracture has an elliptical geometry [10|; a schematic of the system is 
represented in Figure [H 

High pressure in the aquifer pushes the fluid into the fracture, which 
may cause it to further propagate. There are several physical mechanisms 
controlling the behavior of a fluid-driven fracture. For a typical well-driven 
hydraulic fracturing operation, the injected flow rate is high and the fracture 
propagation rate is limited by two dissipative processes: fracturing of the 
rock (controlled by the rock toughness) and pressure dissipation in the fluid 
(controlled by fluid viscosity) [ll|. However, the case considered here differs 
substantially. The source of fluid is the aquifer, which has low permeability, 
and therefore the outflow of fluid from it is slow. The rate of fracture prop- 
agation cannot be faster than the flow of fluid that causes this propagation. 
Since the fluid outflow from the aquifer is the rate-limiting process for the 
fracture propagation, it is the only process considered below. 

If a fluid-driven fracture propagates in a permeable media, the fluid may 
seep into the rock through the walls of the fracture. When permeability of 



the rock is high, this effect may noticeably affect the rate of propagation |12| . 
but in our case a fracture propagates in shale with typical permeabilities of 
the order of 10~^ mD 13|, so the leak-off effects can be neglected. 



We denote the fracture length L and the aperture (maximum width) w, 
which are both functions of time t when the fracture propagates. We assume 
that the initial length of the fracture L{Q) is negligibly small compared to 
its length as it propagates. When describing hydraulic fractures different 
authors have used different relations between w and L. For example, 10| 



and [iJI used w oc L and [15i] and \X^ used w oc y/L. We represent these 
dependences as 

w{t) = aL^{t), (1) 

where a is the constant multiplier and 7 is either 1 or ^ depending on which 
model we consider. The results derived below are applicable to both cases. 

We note by y the direction into the aquifer, so that —y is the direction 
of fracture propagation and |/ = denotes the interface between aquifer and 
caprock. At each time the fracture is filled with fluid, which originated in 
the aquifer and flowed into the fracture. Therefore the length and aperture 
of the fracture satisfy the material balance equation for fluid 

^ = Mt)g(t)|^^„, (2) 

where V{t) is the fracture 'volume' and q{t) is the Darcy flux. Assuming an 
elliptical geometry of the fracture [lC|, we have V{t) = jw(t)L(t). Then, 



using Eq. ([T]), and Darcy's law, we can rewrite Eq. ([2]) as 

dL{t) 4 k dp{y,t) 



dt 7r(7 + 1) /U dy 



(3) 

y=0 



Here k is the aquifer permeability and /i is the fluid viscosity. Therefore, in 
order to determine the evolution of the fracture length L{t) we need to solve 
for the pressure distribution in the system consisting of the aquifer and the 
fracture. In the model problem considered here, the pressure distribution 
can be considered in only one, vertical, direction (see Figured]). 

Below when making approximations and estimates, we use the geome- 
chanical parameters and material properties for the aquifer on the Krechba 
field (Joint Industry Project, In Salah, Algeria) summarized in Tabled] with 
the corresponding references. 



3. Analytical Solution 

The complex multi-physics problem of fracture propagation could be re- 
duced to consideration of the rate-limiting process because of the low perme- 
ability of the aquifer, which allows us also to solve the problem analytically. 
In particular, we assume that the permeability of the aquifer k is much lower 
than the permeability k inside the fracture, i.e. 

k-^k. (4) 

When the fracture propagate its aperture w{t) increases in time according 
to Eq. ([T]). The increase of aperture causes the increase of permeability 



/c = A;(t) ~ w'^{t)/ 12 17|. However it makes the strong inequality @ even 
stronger. For A; = 50 mD the strong inequality (jlj) is valid when the fracture 
aperture w > 2.5 x 10~^ m. Using the parameters from Table [H our estimates 
for the initial fracture aperture based on the work [10] give the initial fracture 
aperture w{0) > 10~^ m, so the strong inequality (jl]) is fulfilled from the 
fracture initiation. 

A fracture propagates if the fluid pressure at its tip exceeds the confining 
total stress |10|. Therefore, we need to consider the pressure change both in 
the aquifer (providing fluid outflow) and in the fracture (providing fracture 
propagation). Below we consider the problem in two approximations. First 
we assume that due to the strong inequality (jlj), the pressure in the fracture 
is established instantaneously, so pressure change is considered just in the 
aquifer. In the second approximation, while still employing strong inequal- 
ity dl]), we do not assume instant equilibration of pressure in the fracture. 
Thus, the pressure equations have to be solved in both the fracture and the 
aquifer, which complicates the problem. Another complication is related to 
the boundary condition at the moving tip of the fracture. However, using 
the self-similarity of the problem, we are still able to obtain an analytical 
solution. 

3.1. Solution with a Constant Pressure at the Fracture Opening 



The pressure evolution is governed by a diffusion equation [18 1 



dt ' dy 

where 



dp{y,t) d'^p{y,t) 






is the diffusion coefficient for pressure, is the porosity of the aquifer, and 
C is the compressibihty of the fluid. 

According to our assumptions the initial and boundary conditions to the 
diffusion problem are the following: initially the pressure is uniform, 

v{y.t)\t=o = vo, (7) 

the pressure in the aquifer far from the fracture remains constant 

v{y.t)\y=oo=vo. (8) 

and the pressure at the opening of the fracture (and inside the fracture) is 
equal to the confining stress in the caprock 

p{y.%^o = ^- (9) 

In such a formulation the diffusion equation is solved only inside the aquifer 
(y > 0), and the role of the fracture is reflected only in the boundary condi- 
tion Eq. (P. 

The solution of this diffusion problem is given by |19i] , 

y 

Calculating the partial derivative of pressure using the solution (11 01) we obtain 



p{y, t) = a + {po- a)erf 



(10) 



dp{y, t) (po - 0-) 



exp 



r 



Acft 



dy 2A/c/t 

Finally, substituting Eq. ( ITTl) into Eq. (jS]), we have 

dL{t) 



/H 



4 






:ii) 



12) 



dt 7r(7 + 1) 

Integration of the latter equation in time, with the initial condition -L(O) = 0, 

m = /3otl/^ (13) 



gives 
where 



/3o 



7r(7 + l)V/^ 



(14) 



When there is no limitation related to the slow outflow, propagation of hy- 
draulic fracture in non-permeable rock has a different exponent: | 20|. The 



square-root time dependence has been predicted for a hydraulic fracture prop 



agation with a large fluid-loss rate [16 



3.2. Coupling the Aquifer and the Fracture 

Once we do not use the assumption of instantaneous pressure equilibra- 
tion inside the fracture, we have to consider its evolution in both regions - 
the aquifer and the fracture. Since the permeabilities of these regions differ 
significantly [Eq. (|4])], the pressure diffusion coefficients differ as well; there- 
fore the diffusion problems should be considered separately. We start from 
the aquifer {y > 0). 

While the diffusion equation (E]) still holds, the boundary conditions are 
subject to change. The pressure value at the aquifer-fracture boundary de- 
viates from the value at the tip, so boundary condition (|9]) does not hold 
anymore. Therefore, we solve Eq. ([5]) with conditions ([7]), ([8]) and the out- 
flow condition provided by material balance, Eq. ([3]), which is equivalent 
to 

(15) 



dp{y, t) 
dy 


^dL{t) 
y^o dt 


7r(7 + l)/i 
~ A k' 



where 

rrf'-v -I- Tl /; 

(16) 

Since boundary condition ( fTSl) includes L(t), it couples the diffusion problem 
in the aquifer with the one in the fracture. Thus, the diffusion problem in the 
aquifer cannot be solved separately, unless additional assumptions are made. 
The results of Subsection 13. II revealed that when the pressure in the fracture 
changes instantaneously, the fracture grows proportional to the square-root 
of time [Eq. (!T3ll]. In the current Subsection we use the same assumption, 

m = fit^i\ (17) 



which will lead us to a self-similar solution of the probleno. Note that the 
coefficient /3, which determines the fracture propagation rate, may differ from 
the coefficient /3o in Eq. [T31 Therefore, calculating the coefficient /3 is our 
main goal. 

The diffusion equation in one space dimension, with constant coefficients, 
is self-similar |2l| , and can be rewritten in terms of one dimensionless variable 
proportional to yj^fcji. Eq. (TT7|) makes the boundary condition ( IT5l) self- 
similar as well, thus the entire problem can be reduced to a single variable. 



^Eq. ()17p and even the expression for parameter /3 (except for a numerical multiplier) 
can be obtained from a scaling analysis of the problem, presented in the Appendix. 



The convenient dimensionless variable is 

P^yim, (18) 

then 

p{y,t)=p{p). (19) 

Using Eqs. (fT7|) - f lT9|) . the diffusion equation ([5]) can be rewritten as 

'iM + f/m . 0. (20) 

Eq. ([8]) gives boundary condition 

P{p)\p=oo=Po (21) 

and, substituting Eq. flT7|) and Eq. flTSj) into Eq. flTSj) . we obtain 

dp{p) 



dp 



= b'- (22) 



Eq. (120]) with boundary conditions ([21]) and ( 12^ gives a self-similar pressure 
profile in the aquifer 

P(P) = Po - ^/3v/^erfc (^p) (?/ > 0), (23) 

where the coefficient /3 is not known yet. 

The pressure diffusion in the fracture is determined by the equation 

dpjy, t) _ d^pjy, t) 

where coefficient Cf differs from the corresponding parameter in Eq. (^. 
Unlike the diffusion coefficient c/, the pressure diffusion coefficient c/ inside 
the fracture changes in time c/ = c/(t) because permeability of the fracture 
k = k{t) increases with the fracture aperture. 

The fracture is in mobile equilibrium, therefore the boundary condition 
at the tip is 

Piy,t)\y=.Li,^ = cJ. (25) 



Another boundary condition, necessary for solving Eq. (1241) is Eq. fllSp . 
Since the assumption (1171) makes the boundary condition fl2^ consistent 
with self- similarity, then analogous to the diffusion problem in the aquifer, 
considered above, the problem flMl) . (1251) and flT5|) can be rewritten in terms 
of the variable p [Eq. (IT8|) ]: 

rfV(p) ^ P' dp{p) 

+ TT^P — -, = 0, 



dp' 



2c f dp 



(26) 

(27) 

and Eq. ([221). Eq. ((261) with the boundary conditions ([27]) and ([22]) leads to 
the following solution 



P{p)\ 



p=-\ 



a 



bd 

P{P) = ^+2^^^ 



erf 



/3 



2^^ 



K,p 1 + erf 



/3 



2^/^ 



(28) 



where the parameter k is defined as 



K = 



The pressure in the aquifer (y > 0), given by Eq. 

y = 0~^ has the value 

b 

P(O^) =P0- ^Py/T^, 



(29) 

at the boundary 

(30) 



while the value given by the pressure in the fracture [Eq. f l25]) ] at y = is 



P(0 ) = a + --y/W^eii ( 



/3 



-K 



2v^ 



(31) 



The continuity of pressure at y = provides p(0 ) = p(0+), which, in accor- 
dance with Eqs. ( l30l) and ( 13T]) . gives 



Po-(r = ^f3y/7icj 



1 + K'^eii 



(3 



-K 



2v^ 



Eq. ( 132]) can be rewritten as 



^ [1 + K-^erf «)] = 5, 



(32) 



(33) 



where ^ is the new dimensionless variable defined as 



e = ,r^. (34) 



and 5 is the constant dimensionless parameter 

^ _ (po - a) 



(35) 



To solve the transcendental equation (!33|) we need to estimate character- 
istic values of parameters k and 5 for our system. Eqs. (135|) . (TT6|) . ([6]) and 
Table[T]give 5 = 10~^ and 5 = 1.5 x 10~^ for 7 = | and 7 = 1 correspondingly, 
therefore 

(5 < 1. (36) 

The pressure diffusion coefficient in the fracture c/ is determined by the 

expression similar to Eq. ([6]), so k = yr^, which provides 

K < 1. (37) 

even when the fracture porosity (j) = 100%. Inequality (!37|) becomes even 
stronger with time due to the increase of fc = k{t). Since the error function 
is positive for any positive argument, the expression in the square brackets 
on the left-hand side of Eq. (l33ll is > 1. Therefore, the strong inequality 
fl36|) provides ^ -C 1, then, due to Eq. fl37|) . k^ <t^ 1 and, thus, erf(K^) ^ 1. 
Finally we have ^ — S, which, using Eqs. flMj) . f p5|) and f lT6|) . gives 




/3 = — ^p^— iJ-C0. (38) 



The latter expression differs from Eq. (fT4l) only by the numerical coefficient. 
Note that since Eq. (!38l) does not contain k, it does not depend on time. 

Figure [2] presents the results of calculations of the fracture length evolu- 
tion L{t) obtained using Eq. ( iT3l) and Eq. (ITTj) with the values of physical 
parameters from Table [1] Note that Eq. f lT3|) predicts slightly slower frac- 
ture propagation than Eq. fITTj) . The assumption of constant pressure in the 
fracture, used when deriving Eq. f lT3|) . leads to lower pressure gradient at 
the aquifer-fracture boundary; to slower outflow; and therefore to a lower 
propagation rate. 
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It is worth noting that the fracture propagation, predicted by Eqs. (IT7I) 
and ( 155]) . does not depend on parameter k, i.e. on the diffusion coefficient Cf 
in the fracture. This is in hne with the initial understanding of the model: 
the fracture propagation is driven and is limited by the flux from the aquifer, 
which has low permeability. If permeability of the aquifer were high, the 
condition fl57|) would be violated, the solution of Eq. fl5^ would depend on 
K, so the resulting expression for (3 would depend on n and therefore on Cf. 
Besides that, if the permeability of the aquifer were high, the problem would 
not be self-similar (see the Appendix). 

3.3. Correction for the Depth 

Although in our schematic in Figure[T]the fracture is vertical, the solution 
we derive is applicable to fracture propagating in any direction. In the current 
subsection we consider the vertical fracture only. 

The driving force for fracture propagation is the difference between the 
fluid pressure in the fracture and the confining total horizontal stress in the 
caprock. In the derivation above we assumed that the confining stress does 
not change with the depth. We also did not include the hydrostatic change in 
fluid pressure. Both effects can be important when the fracture propagates 
high enough toward the earth surface. Let us estimate how these values vary 
with the depth, i.e. consider p = p{y) and an = crniy)- The values of y are 
non-positive for the consideration below (Fig. [T]). 

The total horizontal stress in the caprock (Juiy) is by definition the sum 
of the effective horizontal stress o"^(i/) and the water pressure Pw{y) in the 
caprock 

(yniy) = (^'niy) + Pw{y)- (39) 

Here we use the sign convention for soil mechanics: compressive stresses 
have positive values. The caprock is saturated with brine (water) and CO2 
from the aquifer does not enter it; also the pressure in the caprock is not 
perturbed by the high pressure in the aquifer due to low permeability of 
shale. The effective horizontal stress is related to the effective vertical stress 
(Tyiy) through the lateral stress coefficient K^ 

a'M = KHcr'yiy)- (40) 

The vertical effective stress, in its own turn, can be easily calculated: 

^viy) = <(0) + (1 - <Ps){Ps - P^)gy, (41) 

11 



where cry(O) corresponds to the caprock-aquifer boundary, 0^ and Ps are the 
porosity and density of the caprock respectively, p^ is the density of the brine 
(water), and g is the gravitational acceleration. Evidently, according to Eq. 
PTj) the effective vertical stress decreases with elevation, since the possible 
values of y are non-positive. 

The fluid pressure in the fracture at a certain height y is given by 

Pf{y)=pA^) + pf9y. (42) 

where pj is the fluid density, calculated in accordance with the value of 
residual saturation (Table [1]). The water pressure in the caprock also changes 
with the depth 

Pwiy) = PwiO) + Pw9y- (43) 

Finally, collecting Eqs. fl5U]) - f H5]) . we obtain the dependence of the driving 
force on the depth 

Pfiv) - crniy) = PfiO) - (Th(0) - [fC^(l - 0,)(p, - p^) + {pj - p^)] gy (44) 

where P/(0) — cr//(0) = po — a and the last term is a positive value, growing 
with elevation (due to y < 0). Substituting the values of the physical pa- 



rameters for our system, and using Kh = 0.46 [22|, calculated based on in 



situ stresses reported in |23[, we obtain the increase of the driving force per 



1 meter of decrease of the depth 

A[pf{y)-aH{y)] = 7.1kFa/m (45) 

Introducing depth-corrections to the initial formulation of the problem 
would change its dimensional picture, so that there will be no self-similarity. 
However, we can take it into account by including corrections to the driving 
force in the resulting expressions obtained using the self-similar solution. 

4. Conclusion 

The safety of CO2 storage in deep saline aquifers relies on the integrity of 
the caprock; fractures in the caprock may serve as pathways for CO2 leakage 
if they propagate long enough. In this paper we present a theoretical model 
for propagation of a fracture driven by fluid outflow from a low-permeability 
aquifer. Due to the slowness of flow, we are able to reduce the problem to 
the solution of the pressure diffusion equation in a system consisting of the 
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aquifer and the fracture. While the problem is complicated by the boundary 
condition on the moving boundary at the tip of the fracture, we are able to 
find an analytical self-similar solution. Based on our solution of the diffu- 
sion problem, we derive an analytical expression for the fracture propagation 
length as a function of time. Our simple model can be used together with 
the results of geomechanical simulations when the explicit consideration of 
fracture propagation is not included. This approach provides an estimate to 
the rate of fracture propagation based on the results of continuum mechanics 
simulations, without involving laborious simulations of fracture propagation. 

Using the geomechanical and material parameters for the aquifer at In 
Salah, we predict the length of a hypothetical fracture propagation to be of 
the order of 25 — 35 m within 10 years after initiation. We also estimate 
the depth correction to the driving force for the fracture propagation. We 
show that the changes of confining horizontal stress and hydrostatic pressure 
with elevation lead to additional increase of the driving force for fracture 
propagation of the order of 7 kPa per each meter of elevation. Besides the 
In Salah site, the proposed model is also applicable to a number of aquifers 
currently used for CO2 storage. As such, the following sites have aquifers with 
low permeability (~ 1 — 10 mD): Nagaoka (Japan), Alberta Basin (Canada), 
MRCSP Michigan Basin (USA), Gorgon (Australia) 0- 

Although our estimate predicts relatively slow propagation, fracturing of 
the caprock can still be a serious safety concern. In order to arrest the frac- 
ture, the fluid pressure must be decreased, which places serious constraints 
for safe CO2 injection in aquifers. 
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Appendix: A Scaling Analysis of the Problem 

Analysis of the main dimensional parameters, which govern the problem, 
allows us to obtain a square-root dependence of fracture length on time, Eqs. 
( 1T7|) and (155]1 . except for a numerical multiplier. Let us rewrite the diffusion 
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problems in dimensionless form, introducing dimensionless pressure as 

fM^^^^^l^. (46) 

Po-cr 

Eq. (|5]) gives 

^-,4^, (47) 

Eq. ([24D gives 

while boundary conditions (|H]), ^^Si), and (TTHj) are equivalent to 

fiy,%^oo = h (49) 

fiy,%=o = o (50) 

and 



5/(1/, t) 



9|/ 
where 



- ^.^. (81) 



P^ ^ 4 1 *^_ 

7r7 + l/i Po-cr 

is a constant, which has the dimensions of a diffusion coefficient. 

Since, according to the strong inequality (^, diffusion in the fracture, 
determined by Eq. ( HHl) . is much faster than diffusion in the aquifer, the 
solution in the fracture is slaved to the time-dependent value in the aquifer. 
Therefore, we can relate dimensional y and t by the (slow) diffusion Eq. ( H7|) 
with the corresponding diffusion coefficient c/ 

y oc {cfty/'. (53) 

The fracture length is determined by Eq. (!5T|) . which suggests a scaling 
relation 

L{t) oc ^. (54) 

Substituting the scaling estimate Eq. ( l53l) into Eq. ( 15^ . we obtain 

m ex ^tv^ (55) 

which is equivalent to Eqs. ( 1T7|) and (l38l) . except for the numerical prefactor. 
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Table 1: Properties of the site for CO2 injection 

Permeability of aquifer* k 5 x 10"^^ m^ 

Permeability of caprock^ ks 10^^^ — 10^^^ 

Porosity of aquifer^ 0.17 

Porosity of caprock^" 0^ 0.01 

aquifer temperature* T 90°C 

CO2 injection temperature* Tj 50°C 

Density of CO* pcoa 900 kg/m^ 

Density of brine^ p^ 1000 kg/m^ 

Density of caprock* ps 2400 kg/m^ 

Compressibility of CO* 1.3 x lO^^ Pa 

Compressibility of brine^ 4.1 x 10^^*^ Pa^^ 

Viscosity of CO^ 9.0 x 10"^ Pa-s 

Viscosity of brine^ 3.0 x 10~^ Pa-s 

Fluid pressure in aquifer* po 30 MPa 

Confining horizontal stress* a 28 MPa 

Water residual saturation* 5*,.^ 0.25 

* from S 

^ from 9[| and [5|] 

§ from 23], values for water are used 
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Figure 1: System of an aquifer and caprock with a vertical elliptic fracture in it. The fluid 
flux from the overpressured aquifer drives the fracture propagation. 
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Figure 2: Predictions for the fracture length evolution, calculated using Eq. (|13p and 
Eq. (J17p for possible values of parameter 7 and physical parameters from Table [1] The 
parameters of the mixture were calculated using the average weighted with the saturations. 
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